分子動力学シミュレーションコードISISの開発

神林 奨

JAERI-M 92-080, 45 Pages, 1992/06

JAERI-M-92-080.pdf:1.3MB

分子動力学法を用いたコンピュータシミュレーションは、物理、化学、生物学など広範な分野において、実験測定や理論計算と相補的な研究手法として定着しつつある。分子動力学法は、ニュートンの運動方程式を数値的に解くことにより、多数の古典粒子からなる体系の運動の様子を数値的に再現する手法である。1980年代に入り、分子動力学法の技術的改良がなされ、実験測定との比較に便利な温度や圧力が一定という条件でのシミュレーションが可能となった。本報告では、まず、分子動力学法の概要をまとめ、次にソフト球モデルを対象として新規に開発した分子動力学シミュレーションコードISIS(isokinetic simulation of soft-spheres)について述べる。ISISコードは、一般的なミクロカノニカルアンサンブルのシミュレーションとともに、束縛法を用いた温度一定のシミュレーションを可能としている。また、運動方程式の逐次積分に第5次の差分式を用いて数値誤差を小さく抑え、粒子間力の計算時間を短縮するために簿記法を用いている。